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SUMMARY 


A general Rayleigh analysis is used as a basis for developing four 
methods of flutter analysis that are applied to twelve low-aspect-ratio 
wings. These wings were previously tested at a Mach number of 1.3 Ly 
progressively varying certain wing parameters until flutter occurred. 

They were rectangular in plan form and had aspect ratios between 5*00 
and 4 . 55 . The four methods of flutter analysis used are: section coef- 

ficients for harmonically pitching and translating rectangular wings in 
a Raylei^ type of analysis, two-dimensional coefficients in a Rayleigh 
type of analysis, total coefficients for harmonically pitching and trans- 
lating rectangular wings in a representative -section analysis, and two- 
dimensional coefficients in a representative -section analysis. Each of 
the four methods involved two degrees of freedom, namely, first bending 
and first torsion of a cantilever wing. 

The analytical results are compared with the previously obtained 
experimental values. The comparison Indicates that the use of section 
aerodynamic coefficients derived on the basis of three-dimensional flow 
leads to a significant improvement in the correlation of theory and 
experiment. 


INTRODUCTION 


The problem of theoretically determining the flutter characteristics 
of unswept wings of low aspect ratio in svg)ersonic flow has become of 
increased interest. Most of the previous analytical work on this problem 
has been based on air -force and moment coefficients for two-dimensional 
supersonic flow, such as those tabulated in reference 1. For example, 
reference 2 presents the results obtained at a Mach number of 1.5^ Ly 
using two-dimensional coefficients in a representative-section type of 
flutter analysis, for twelve \mswept wings with aspect ratios ranging 



2 


NACA TN 5501 


from 5*00 to As explained in reference 2, these wings were also 

tested at a Mach number of 1.3 hy progressively shifting their centers 
of gravity and elastic axes and modifying their bending and torsional 
frequencies until flutter occurred. A comparison of the calculated and 
experimental results showed that in the majority of cases the calculated 
flutter speeds were considerably below the experimental flutter speed. 

This discrepancy suggests in part that, at least in the low supersonic 
speed range, two-dimensional coefficients are inadequate and more real- 
istic aerodynamic coefficients should be used in the flutter analysis 
of unswept low-aspect-ratio wings. 

In reference 3^ streamwise section and total air-force and moment 
coefficients expanded to the seventh power of the frequency of oscilla- 
tion were developed for harmonically pitching and translating rectangu- 
lar wings moving at supersonic speed. The section coefficients were 
used in a Rayleigh type of flutter analysis to calculate the flutter 
speeds of a rectangular wing of aspect ratio ^»55 at several Mach num- 
bers in the low supersonic speed range. For comparison the wing was 
also analyzed by using the two-dimensional coefficients of reference 1 
in a Rayleigh type of analysis. Examination of the results showed the 
flutter speeds based on the rectangular-wing section coefficients to be 
higher than those based on two -dimens ionail coefficients, particularly 
at the lower Mach numbers. Application of a Rayleigh type of analysis 
involving the section coefficients of reference 3 to the wings of ref- 
erence 2 might therefore be expected to yield a better correlation between 
theory and experiment than was obtained in reference 2. 

Also of interest is reference 4 in which a comparison is made between 
flutter results obtained by using two-dimensional coefficients in a 
representative-section type of analysis and total coefficients for rec- 
tangular wings in the same type of analysis. For wing parameters in the 
range of those given in reference 2, reference k also shows an increase 
in calculated flutter speed resulting from the use of finite-wing 
coefficients. 

In the present paper four methods of analysis are applied to the 
twelve wings of reference 2 and the results are compared with the exper- 
imental results in reference 2. These four methods of flutter analysis 
are: section coefficients for a pitching and translating rectangular wing 

in a Rayleigh type of analysis, two-dimensional coefficients in a Rayleigh 
type of analysis, total coefficients for a pitching and translating rec- 
tangular wing in a representative -section type of analysis, and two- 
dimensional coefficients in a representative-section type of analysis. 
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SYMBOLS 

A aspect ratio, s/b 

b one -half chord 

c speed of sound in undisturbed medium 

8h^Sa first bending and first torsion damping coefficients, 
respectively (see ch, IX of ref. 5) 

h vertical displacement of axis of rotation Xq, positive 

downward 

E generalized coordinate in bending degree of freedom, hoe^^^ 

ho bending amplitude at tip of wing 

k reduced frequency, oib/V 

Ih^i^h coefficients of section lift and moment, respectively, 
associated with mode shape 

coefficients of section lift and moment, respectively, 
associated with mode shape 

Li, Mi coii5)onents of section force and moment coefficients, respec- 

tively, for rectangular wing (see ref. 5) ii^ equation (7) 
and for two-dimensional wing (see ref. l) in equation (8); 
i = 1, 2, 5 j an^i ^ 

Li, Mi components of total force and moment coefficients, respec- 

tively, for rectangular wing (see ref. 3); i = 1^ 2, 3> 
and 4 

M Mach n-umber, V/c 

Mq^ aerodynamic section moment on wing about axis of rotation Xq, 

positive leading edge up 

P aerodynamic section normal force, positive downward 

r<x nondimensional radius of gyration of wing section about 

elastic axis, where is mass moment of 

inertia per unit span about elastic ajcis and m is mass 
of wing per unit span 



4 


MCA TN 5501 


s 

t 

V 

X 


Xa 

Xo 

y 

Zh 

Za 

a 

a 

°o 

P = )/iM^ 

K 

P 

(U 

tjoa 


one -half span of wing 
time 

velocity of flow 

nondimensional chordwise coordinate measured from leading 
edge, referred to wing chord 2b 

location of center of gravity of wing measured from elastic 
axis (see ref. l) 

chordwise position of axis of rotation of wing (elastic axis) 

nondimensional spanwise coordinate measured from midspan of 
wing, referred to wing half-span s 

first bending mode shape of wing 

first torsion mode shape of wing 

angle of attack, positive leading edge up 

J yjY*- 

generalized coordinate in torsion degree of freedom, Ooe-*- 
torsion amplitude at tip of wing 
1 

density par ame ter , rtpb m 
density in undistxirbed medi\im 
frequency of oscillation at flutter 
first bending frequency of wing 

first torsion frequency of wing 


METHODS OF FLUTTER ANALYSIS 
Rayleigh Analysis 


General considerations .- The wings to be analyzed are rectangular 
in plan form and were tested as cantilevers in the Langley supersonic 
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flutter apparatus (a 9-inch by l8-inch supersonic dravdown tunnel). In 
a Rayleigh type of analysis of such wings, the bending component of the 
flutter mode can be approximated by the first bending mode of a \uiiform 
cantilever wing and the torsion component by the first torsion mode. 

The flutter determinant is then formed and is solved for the flutter 
condition. (a detailed discussion of the Rayleigh type of analysis as 
applied to flutter may be found in ch. IX of ref. 5* ) 

The bending component h and the torsion component a of the 
flutter mode may be written as 


h(y,t) = Zh(y)h(t) 
a.(y,t) = Zo,(y)a(t) 


( 1 ) 


where y is the nondlmensional coordinate shown in figure 1, and 

are the first bending and first torsion mode shapes shown in fig- 
ure 2, and E and a are the generalized coordinates in the bending 
and torsion degrees of freedom^ respectively. The section aerodynamic 
force or aerodynamic force per unit span, positive downward, associated 
with equations (l) may be written as 


P = 


-4pb2o)2 |^(y)h + 2a(y)^ 


( 2 ) 


and the section moment, positive leading edge up, about the arbitrary 
axis of rotation x = Xq may be written as 


Mq, = -^pb^oj^ I%(y)h + I%(y)a 


( 5 ) 


where o) is the frequency of oscillation, b is the one-half chord of 
the wing, and i% are complex coefficients of the lift and moment 

associated with the mode Z^, and Zq, and are the complex coeffi- 

cients of the lift and moment associated with the mode Z^. Each of the 
aerodynamic coefficients Za^ and m<x^ in addition to being a 

function of the spanwise variable y, is a function of Mach number M 
and reduced frequency k = bo)/V. Although these coefficients may be 
taken to apply at either subsonic or supersonic speed, the present paper 
is concerned only with the supersonic speed range. 
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The equilibrium equations at flutter may be obtained by setting up 
the potential and kinetic energies and the work of the applied forces, 
introducing equations (l), (2), and (3) and the mass and stiffness prop- 
erties of the wing, and then applying Lagrange's dynamical equation, as 
shown in chapter IX of reference 5* From the equilibrium equations a 
flutter determineint may be obtained in the form 


^hh 

^ha 


Aooa 


= 0 


(^) 


where the determinant elements are 


Ahh = 


(^) fo 


k 

^ha “ ^ / ^h^ dy - — K / Zo,Zh dy 

lo ^ Jo 


4 

Aah = Xa Z^Zh ~ ~ 


Vi = V 


1 - f— Y 

\ 05/ 


r.l 


Za^dy 


" Jo 


maZa 


(5) 


From equations (4) and ( 5 ) fo\ir methods of analysis are obtained by using 
various approximations in evaluating the integrals of equation ( 5 )* 


Section coefficients for rectangular wing ,- The following approxi- 
mate expressions for the section coefficients are employed: 
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where Lj^ and (i = 1, 2, 5^ sind are the components of the sec- 

tion coefficients given in reference 5 i’oi* sl rectangular wing in super- 
sonic flow oscillating harmonically as a rigid body in pitch and verti- 
cal translation, (a preliminary unpublished analysis, based on parabolic 
bending of a rectangular wing which closely resembles the mode shape 

suggests that the results obtained by using the distributions of lift and 
moment for the mode shapes and would be nearly identical to the 

results obtained by using the approximate distributions given by eqs. (6), 
when multiplied by the mode shape Z^ or Z^ and integrated in the man- 
ner required in eqs. (5)*) Upon substituting equations (6) into equa- 
tions (5), the determinant elements of equation become 
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The uncoupled first tending mode shape and the first torsion mode 
shape Z^ needed for the evaluation of the integrals of equations (j) 
are shown in figure 2 . The integrals of equations (j) containing only 
mode shapes can be evaluated to give 

Zt^dy = 0.25 




0.557 



0.50 


A numerical method for evaluating the integrals of equation (7) involving 
the aerodynamic coefficients and is given in appendix B of ref- 

erence 3. (in \islng ref. 5 , note that the spanwise coordinate y of the 
present paper and the spanwise coordinate ^ of the reference paper are 
related hy y = 1 - ) 


Coefficients for two-dimensional wing .- If two-dimensional air-force 
and moment coefficients are used in place of the section coefficients of 
reference 3^ the force and moment coefficients in equations (7) appear 
as constants in the integrals and can be factored from under the integral 
signs, and the determinant elements of equation (k) become 
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where and Mj^ (i = 1, 2, 3^ Q-nd 4) now refer to components of 

two-dimensional coefficients, such as those tabulated in reference 1. 


Representative-Section Analysis 

Total coefficients for rectangular wing * - By applying mean-value 
theory to the integrals in equation ( 7 ) and, in the process, by assuming 
the representative section to be the same for all integrals involved, 
the determinant elements can be written as 



where the subscript r denotes evaluation at a representative spanwise 
station y = r. Since the q^uantities having the subscript r cancel 
in the solution of equation (4), equations ( 9 ) may be rewritten as 

Ahh - 1 - - I «(£i + lia) 

Aah = ^ + IW2) 

Aaa = 1 - ^ 
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where L = T dy and ^ = f dy (i = 1, 2, 3, and ^4-) are the 

Jo Jo 

components of the total force and moment coefficients for a pitching and 
translating rectangular wing given in reference 

Coefficients for two-dimensional wing .- If infinite aspect ratio is 
substituted into the aerodynamic coefficients of equation (lO) (see 
ref. 3), the determinant elements can be written as 



where, as in equation (8), Lj^ and refer to components of two- 

dimensional coefficients, such as those tabulated in reference 1. 


Solution of Flutter Determinant 

The flutter condition is determined from the nontrivial solution 
of equation (ij-) obtained by using as determinant elements the various 
approximate forms of equations (5) given by equations ( 7 ), (8), (lO) , 
and ( 11 ) • This condition, which requires that the real and imaginary 
parts of equation (4-) vanish simultaneously for the same set of aero- 
dynamic and wing parameters, may be obtained by various means (see 
ch. XIII of ref. 5 ). 

In the present paper the ratio in equation (4) is replaced 

by the equivalent quantity (oh,/ta)* Then, for a particular wing 

and Mach number, for which values of M, k, Xq, Xq,, and 
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are specified, equation (4) contains the two unknown parameters 

and k = bo)/v. The reduced frequency k (upon which the various aero- 
dynamic terms are dependent) is varied until the same value of Cjo^/o) 

is obtained from both the real and imaginary parts of equation (4). This 
is the required condition and yields the values of k and at 

flutter and consequently the flutter-speed coefficients v/bo^ for the 
wing at the selected value of M* 


APPLICATION AND DISCUSSION OF RESULTS 


The four methods of analysis outlined in the previous section were 
applied to the twelve wings of reference 2. The wing parameters needed 
in these analyses and a description of each wing profile, obtained from 
reference 2, are given in table I of the present paper. The flutter 
parameters, and and consequently calculated by 

these methods are listed in table II. For comparison table II also 
includes the experimentally determined flutter parameters given in 
reference 2. 

In figure 5 the data of table II are plotted in line-graph form. 

The line-graph method of plotting is en 5 )loyed to achieve a separation 
of the data and ease of comparison not otherwise obtained because of 
the insufficient range of variation of the different wing parameters. 

Also shown in figure 5; as flagged points, are the analytical results 
of reference 2. These results were obtained by the last method of the 
previous section (two-dimensional coefficients in a representative- 
section analysis) but included structural damping. Structural dairp)ing 
covild also have been included in the calculations of the present paper by 

replacing by and (X3^ by + iga)> where gj^ 

is the damping coefficient in bending and g<x is the damping coefficient 

in torsion, in the methods discussed previously. Since damping was not 
included, the calculations of reference 2 may serve to indicate the 
effect the inclusion of damping would have on the calculations of the 
present paper. 

Figure 3 (a) shows a comparison for each wing of the values of 
reduced flutter speed V/bco (reciprocal of reduced frequency k) , 
obtained by the four methods of analysis. The results obtained by 
using finite-wing section coefficients in a Rayleigh analysis are 
closest to experiment in all twelve of the cases treated. 

In figure 5(b) values of the ratio of flutter frequency to tor- 
sional frequency o)/o:tc compared for the various wing models. As 
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may be noted in the fig\ire, for five of the wing models (A-1, B-1, C-1, 
C-2, and F-l) comparatively large differences between the theoretical 
and experimental values of exist; these differences woiad probably 

be reduced by the inclusion of more degrees of freedom in the various 
analyses. However, it may be seen by con 5 )aring the flagged and unflagged 
right triangles in figure 3(b) that the inclusion of structural damping 
in the Rayleigh analysis involving section coefficients for a rectangu- 
lar wing may sufficiently reduce the differences between experiment and 
theory . 

Figure 5(c) shows a con 5 )arison of the values of flutter-speed coef- 
ficient V/ba^ determined for the various wing models from the data 

presented in figures 5(a) and 5(b). The resialts of using the Rayleigh 
analysis involving section coefficients for rectangular wings are 
closest to experiment in the majority of the cases treated, that is, 
except for models A-1, C-1, and C-2. The section-coefficient results 
in these cases are above the experimental values (nonconservative). 
Inclusion of more modes in the analysis may relieve this situation. 

Also of interest in the present comparison are the c\arves of ''I^bix^ 
calcinated in reference 5 f’oi’ model B-1 of table I in the Mach number 
range 10/9 < M < 10/6 by the first two methods of the previous section, 

that is, section coefficients for a nondeformable rectangular wing in a 
Rayleigh type of analysis and two-dimensional coefficients in the same 
type of analysis. These curves, taken from figure 12 of reference 5; 
are shown in figure 4. The main feature of these curves, as pointed 
out in the reference paper, is that the use of finite-wing coefficients 
is very influential at Mach numbers near unity but, as would be expected, 
becomes less so as the Mach number is increased. At M = 10/9, for the 
particular wing analyzed, the flutter speed obtained by using two- 
dimensional coefficients is about 62 percent of that obtained by using 
rectangular-wing section coefficients, whereas at M = 10/6 it is about 
95 percent. For comparison at M = 1-5 the experimental value for 
model B-1 and the results of using two-dimensional coefficients with 
and without structural dandling and total rectangular-wing coefficients 
in a representative -sect ion analysis are included in figure 4. The 
values plotted in figure 4 at M = 1.5 are, of course, also given for 
model B-1 in figure 5(c). As may be noted in figure 4, the result 
obtained by using rectangular -wing section coefficients in a Rayleigh 
analysis is in excellent agreement with experiment. 


CONCLUDING REMARKS 


The results of applying four methods of flutter analysis to a 
series of twelve wings have been presented and discussed. The wings in 
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question, which were fluttered previously at a Mach mamber of I .5 in the 
Langley supersonic flutter apparatus, had aspect ratios ranging from 5 . 00 
to 4.55 and various profile shapes, masses, and stiffness properties. 

The four methods of analysis, which are derivable from a general Rayleigh 
type of analysis, are: section coefficients for a pitching and trans- 

lating wing in a Rayleigh type of analysis, two-dimensional coefficients 
in a Rayleigh type of analysis, total coefficients for a pitching and 
translating rectangular wing in a representative -section analysis, and 
two-dimensional coefficients in a representative -section analysis. Each 
of the four analyses involved two degrees of freedom, namely, first 
bending and first torsion of a cantilever wing. The section and total 
aerodynamic coefficients for rectangular wings that were used are those 
that were developed, for wing pitching and vertical translation, to the 
seventh power of the frequency in NACA TN 50T6. 

Previous analyses of the flutter of unswept wings of low aspect 
ratio in supersonic flow have customarily involved the use of aerodynamic 
coefficients for two-dimensional flow. The present paper shows that the 
use of aerodynamic coefficients for rigid-body motions of a wing, namely 
pitching and vertical translation, derived on the basis of three- 
dimensional flow leads, at least in the low supersonic speed range, to 
a significant improvement in the correlation of theory and experiment. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. , August 15; 1954. 
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COMPARISON OF CALCULATED AND EXPEBDffiMTAL FLUTTER PARAMETERS 


(a) v/bo) 


Values of V/bo) 


Rayleigh analysis 


Representative -sect Ion 
analysis 




Rectangular- 
wing section 
coefficients 

Two-dimensional 

coefficients 

Rectangular- 
vlng total 
coefficients 

Two-dimensional 

coefficients 

A-1 

10.15 

8.96 

4.65 

7.09 

4.81 

B -1 

9.98 

9.10 

6.58 

7.89 

6.98 

B -2 

10.51 

8.45 

7.02 

7.39 

7.17 

B -3 

10.20 

8.87 

7.52 

7.78 

7.49 

B-U 

10,40 

9.10 

7.46 

8.10 

7.45 

B -5 

10.03 

9.60 

7.97 

8.25 

7.85 

C -1 

, 9.715 

8.78 

5.45 

7.16 

5.45 

C -2 

i 9.92 

8.90 

5.90 

7.29 

5.93 

D-1 

9.04 

7.10 

4.55 

5- 97 

4.50 

E -1 

19.15 

18.27 

15.48 

16.05 

14.32 

F-1 

19.61 

14.45 

8.12 

11.65 

8.55 

G -1 

7 . 71 

6.95 

4.65 

5.19 

4.64 

(b) u)/t% 


Values of u)/< 


Rayleigh analysis 


Representative-section 

analysis 


Rectangular- 
ving section 
coefficients 


Tvo -dimens ional 
coefficients 


Re ctangul ar- 
ming total 
coefficients 


Two-dimensional 

coefficients 


Values of v/bo^ 


Experiment 


Rayleigh analysis 


Rectangular - 
wing section 
coefficients 


Representative-section 

analyslB 


Two-dimena ional 
coefficients 


Rectangular- 
wlng total 
coefficients 


Two-dimensional 

coefficients 
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Figure 5*- Comparison of calculated and experimental flutter character- 
istics for twelve wing models described in table I. 



Model designation 

B-2 B-3 B-4 B-5 C-l C-2 
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(b) Ratio of flutter frequency to torsional frequency, 
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O Experiment 

Reef, wing section coeff. I Rayleigh 

Two-dimensionol coeff. j analysis 

A Rect. wing total coeff. I Repre s en ta tive 
^ Two-dimensional coeffj section malysis 

A Calc, results of ref. 2, g^ and / 0 


Figure 4.- Flutter-speed coefficients plotted against Mach number for 
model B -1 of table I. {k = 4.53; ^ = 95-5; Xq = 0.54l; = O.35O; 

= 0.59; and = O. 583 ) 
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